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1 Introduction 



In spite of the photon being a fundamental gauge particle of electromagnetic 
interactions, it is known to develop a hadronic structure in its interactions 
with matter. The notion of the hadronic structure function of the photon, 
F^, is introduced in analogy to the well known nucleon case. The first mea- 
surements of F2 became available from e^e~ collisions in which one of the 
leptons is scattered under a small angle. These interactions may then be 
interpreted as processes in which a highly virtual photon, of virtuality Q^, 
probes an almost real target photon, with virtuality ~ 0. 

While the proton structure function has been well measured over a 
wide range of and the Bjorken variable x [1] , data cover a restricted 
kinematic range (0.001 < x < 0.9) and are subject to much larger systematic 
uncertainties [2] . This is mainly due to experimental limitations in measuring 
the centre of mass energy W of the 7*7 system, in particular at large values 
oiW. 

In the present paper, a new parametrization of the parton distributions in 

the photon is extracted in next-to-lcading order (NLO) of perturbative QCD. 
It differs from other NLO parametrizations [3-7] in that the data used in the 
fitting procedure include the expected behaviour of at low-x [8] , as derived 
from F2 measurements [9] under Gribov factorization assumption [10] and, in 
addition, the measurements of the dijet photoproduction cross sections [11]. 

2 Theoretical background 

2.1 Physical quantities and parton distributions 

In QCD, the hadronic cross sections are given as convolutions of coefficient 
functions and parton densities. In higher orders of perturbative expan- 
sion, the definition of parton densities is not unique and depends on the 
adopted factorization scheme. In this paper we adopt the DIS7 factoriza- 
tion scheme [3], where the relation between and the parton distributions, 
f].{x, Q"^) {k = q or G), is given by 
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where F2{x,Q^) = ^F2{x,Q^), a^iQ^) is the QCD running couphng con- 
stant, (g) denotes convolution, A^f is the number of flavours and is the 
electric charge of quark q in units of e. The parton distributions are convo- 
luted with known coefficient functions^ Cf,2{x), Cg, 2(yX) [3,13]. In Eq. ((H), 
all the iVf flavours are assumed to be massless and the effects of heavy quarks 
are discussed in Sec. 12.31 

The cross-section for dijets photoproduction depends on both the photon 
and the proton parton distributions and we take the latter to be fixed by the 
deep inelastic ep scattering measurements. To calculate this cross-section at 
NLO, the program by Frixione and Ridolfi [14] is used. 



2.2 Evolution equations for the photon 

In order to model the parton content of the photon, f^{x, Q"^) {k = q, q, G), 
we parametrize the parton distributions (PDF) at some scale Qo, f^{x,Ql), 
and evolve them to other scales through the analogue of the DGLAP evolu- 
tion equations for the photon which contain an inhomogeneous term [15-17]. 
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where PkA {x, Q^) are the splitting functions, and is the fine-structure 
constant. 

Within next-to-leading order accuracy the splitting functions are given 

by 

PkA (x, Q') = (x) + ^^pS (x) , (3) 

and the higher order terms P^^^(x) depend on the adopted factorization 
scheme. The QCD splitting functions for A = q,q,G at NLO in the MS 
scheme are given in [18]. The photon splitting functions, P^^ {XyQ"^) can be 
obtained from the gluonic ones {Pqc, Pgg) by setting Co = and multiplying 
by factor three. For the DIS7 factorization scheme the latter still have to be 
transformed as described in [3]. 



2.3 Heavy quarks 

In order to solve the evolution equations given by Q we still need to specify 
the treatment of the heavy fiavour contribution. 

^Cq \ used here is the one given in Appendix I of [13] divided by A^f. 
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The approaches traditionally used are either the Fixed Flavour Number 
scheme (FFNS) or the Zero Mass Variable Flavour Number scheme (ZM- 
VFNS). While FFNS is known to be adequate at low values, the ZM- 
VFNS is the appropriate scheme for calculations at high . Recent devel- 
opments [19-22] aim at constructing schemes which work well in the whole 

range. Any such scheme has to approach FFNS at low and ZM-VFNS 
at high (5^. In the following we propose a phenomenological approach which 
smoothly interpolates between FFNS and ZM-VFNS results. 

In general one can write 
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where the contributions from the light and heavy quarks, e^HhiQ^), 
given separately. 

In FFNS there are no heavy quarks in the probed target and a pair of 
heavy quarks can only be produced in the final state. The threshold condition 
for such production is for the hadronic state mass, W, to be above 2mft, i.e. 
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where nih is the heavy quark mass. 'Hh{Q'^) is calculated at fixed order in 
ttem and as{Q'^) with ruh 7^ 0. Within the accuracy needed for the current 
calculation only 7*7 — > hh and 7*G — > hh processes contribute, i.e. Hh — 
Hi + H^. The point-like photon contribution from 7*7 scattering, T-Q, is 



nl{x,Q') = 6el^a,,{x,Q') 



while the contribution from 7*6* scattering, is 
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o-fih is the Bethe-Heitler type reduced cross section [3,4,23], 
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For ^ we can neglect terms vanishing in the hmit m\/Q'^ 
which results in the approximation 
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The large logarithm, In-^, spoils the accuracy of this fixed order calcula- 
tion and in order to achieve the NLO accuracy this term would have to be 
resummed. This resummation would generate the f1 distribution. 

In the ZM-VFNS the heavy quark masses, m/^, serve only as transition 
scales. When crosses the value of m|, the number of active fiavours, A^f, 
changes by one. The effect on the evolution equations is twofold. First, the 
limit of summation over quarks changes and, second, the dependence of 
otsiQ"^) changes. As we want as(Q^) to be continuous across the transition 
scales, we adjust the scale parameter of QCD, AqcolA^f), so that 



a,{ml- Ni, AQCD(A^f)) = a^iml Ni + 1, AqcoliVf + 1)) • 



(10) 



Otherwise there is no dependence on quark masses in the evolution. The 
evolution is done according to Eq. (j21) and Hh is given by the same formula 
as for light quarks. 
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This ZM-VFNS formula sums up large collinear logarithms which were dis- 
cussed in the context of validity of Q- The remaining terms of Q, con- 
tributing to the 7*7 process, are included in the definition in the DIS7 
scheme. The same terms in the 7*G contribution lead to the second term of 
TC^ (see Eq. (fTT|l ). However, at finite the dependence on starts to be 
important and in a more precise calculation [20] it is taken into account by 
the mass dependent coefficient functions. 

In summary, FFNS has no resummation of the large collinear logarithms, 
important at high Q^, and ZM-VFNS has no mass dependence, important 
at low and near the W = 2mh threshold. Note that both ZM-VFNS and 
FFNS schemes are formally correct pQCD results. At intermediate scales 
the contributions from both powers of m\/Q'^ (Eq. (jSJ) and collinear logs 
resummation (Eq. should be taken into account. What we are looking 

for now is a prescription which would be a good approximation over the 
whole range of Q^. 

It is clear from the above discussion that simply adding the two con- 
tributions would result in some double counting. A consistent solution is 
not uniquely given by QCD [24] and several approaches have been dis- 
cussed [7,19,21,22]. Here, we propose to construct the heavy quark con- 
tribution to F2 as a weighted sum of FFNS and ZM-VFNS expressions with 
a Q^-dependent weight, S^v, 

=[l- S^ml Q')] {HliQ') + Hf{Q^)) + S,.{ml Q')ni\Q') . 

(12) 

Intuitively, S^v quantifies the amount of evolution and it will be defined so 
that 5ev = at < m| and iSev — ^ 1 for Q'^ ^ m^. By construction, 
the coefficients of (Hi + H^) and in (fT^ add up to 1 to prevent dou- 
ble counting. Using pi|l we can explicitly write out the quark and gluon 
contributions 

I Zn 'J 

+ (l-5e.)^a,,®/^+5ev^Cg)2®/G- (13) 



The NLO corrections coming from the '~f*G hh process are given by the 
last two terms where ath plays the role of the mass dependent coefficient 
function. As this is an NLO contribution and we have no mass dependence 
in the other coefficient functions, we take 



as 



Hh = (1 - 5ev)H^ + 25ev 1 + ^4^2 ® 

ZTT 



+ ^cSl^fa- (14) 



5 



This last approximation has the advantage of noticeably speeding up the 
minimization procedure, at a cost of disregarding a small correction, well 
within the uncertainty of the NLO approximation. 
We define »Spv as 



for <ml, 
min (l, T{ml, Q"^)) for Q'^ >ml, 



where 



In^ 



(15) 
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is the QCD evolution scale. 

To correct the behaviour of massless coefficient functions close to the W 
threshold, we introduce a phase space suppression factor, (Sps. (Sps must be 
zero at the threshold, W — 2mh, and must go to one for W » 2mh. Sps 
is taken as the ratio of the leading twist T-Q over its asymptotic (Q^ » ml) 
value. 
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The final result for heavy quark contribution to F2 is 
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In the NLO calculation of the dijet photoproduction cross sections, as 
implemented in the Prixione and Ridolfi program [14], the partonic cross 
sections assume massless quarks. The factorization scale is taken as the 
average of the jets. The value of £'t for the data considered in this paper 
is well above rrih and below m^. Hence we expect the heavy quarks mass 
effects to be neghgible. 



2.4 Constraints on the gluons 

The gluon content of the photon can be extracted only indirectly, in particular 
from the evolution of F^. The case of the photon is even more challenging 
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than that of the proton. The contribution of the point-hke sphtting of the 
photon into a qq is at the origin of strong scahng violations, positive for 
all values of x. This decreases the sensitivity to scaling violations due to 
gluon radiation. The difficulty is exacerbated further by the lack of the 
momentum sum rule and of small x measurements where the scaling violation 
is dominated by gluons. 

In the present extraction, the lacking small x measurements of are 
replaced by the low x measurements of the proton F2 transformed as sug- 
gested in [8] . The transformation is based on Gribov factorization [10] which 
relates the total 77 cross section to those of 7p and pp. For low x one can 
thus obtain 

F,-(x,g2) = Ff(x,g2)^J^. (19) 



a, 



pp 



Using the Donnachie and Landshoff [12] (DL) parametrization of the cross 
sections, which gives a good representation of the data, one obtains at large 
W 

F^/a,^ = 0A3Fl (20) 

In the following, we assume that Gribov factorization holds and only assign 
an error to the transformation (j^IH) which arises from the DL parametrization. 
More specifically, we refitted the pp, pp and •yp total cross section data at 
high energy keeping the Regge intercepts obtained by DL. Depending on 
the energy range fitted, the uncertainty on the ratio CT-yp{W) / app{W) varied 
anywhere between 1% and 6% and we chose to use 3% as an error estimate. 

Frankfurt and Gurvich derived an equivalent of the momentum sum rule, 
the FG sum rule [26]^. It relates the momentum fraction carried by partons 
in the resolved photon to the physical e~^e^ "hadrons" cross section. The 
authors estimate that for < 10 GeV^ 



— [dxx [S'(x, Q^, P2 = 0) + G{x, Q^, P2 = 0)1 ^ 1 + — In ■ 

^em J J Svr . 



(21) 

where 

E'(x,g^O)= ^ [f^ix,Q') + f^ix,Q')] (22) 

q=u,d,s 

and G{x,Q^,P^ = 0) = /^(a;,^^). The quark distributions in the photon 
used in [26] are normalized such that S'(x, Q^,Q^) is given by the box di- 
agram contribution. At NLO this normalization corresponds to the DIS7 



similar sum rule in LO was derived by [25]. 
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factorisation scheme, where the subleading box diagram contribution is ab- 
sorbed into the quark distributions. 

This sum rule was not directly incorporated in the fits, however it is used 
to assess the quality of the fitted parametrization. 

In an attempt to constrain further the gluon distribution in the photon, 
the cross section measurements of dijet photoproduction were added in the 
fitting procedure. However we found out that these measurements are dom- 
inated by the contribution of gluons in the proton, while the region most 
sensitive to gluons in the photon is suppressed by kinematic constraints. To 
establish this, we calculated the resolved photon contributions to the total 
dijet cross section from various partonic subprocesses as a function of the 
observed photon momentum fraction participating in the hard interaction, 
^obs_ rpj^g cross sections were calculated using the Frixione-Ridolfi code with 
CTEQ5M proton PDFs and our photon PDFs. The results are shown in 
Fig. E Note that the dominant contribution depends on the quark content 




-obs yObs yObs yObs 



Figure 1: Typical contributions of different parton types, as denoted in the 
figure, to da/dx*^^ for E'x ^ [14, 17] GeV. Hadronization corrections are not 
included. 

of the photon which is fairly well established by the data. We checked 
that this conclusion is not affected by the choice of the PDFs used in the 
calculation. 



3 Parametrization 

In practice we solve the NLO evolution equations numerically in the x space. 
We use a programme [27] which performs the evolution over a two-dimensional 
grid in x and In(Q^). The actual (x, Q^) values are obtained from the grid 
by a quadratic interpolation. 

The dijet cross sections, when included in the fit, are recalculated for each 
set of f2 {x, Q^) and the hadronic corrections are added. 
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Our parametrization of the initial parton distributions, defined at — 
2 GeV^, aims at describing the experimental data below the charm threshold. 
Thus we explicitely parametrize only the u, d, s quarks and the gluon. The 
c, b and t quarks are generated radiatively once their respective thresholds 
(transition scales) are crossed. 

All quark distributions in the photon are parametrized as a sum of point- 
like and hadron-like contributions, 

U^) = f,(.) = <A- \f^,l,\-f!^^ + /r^C-) . (23) 

Apart from the factor, the point-like contribution is the same for all 
quarks. The hadron-like contribution is assumed to depend on the quark 
mass only. For u and d quarks we parametrize it as 

/r°(^) = /r'^i^) = A^^^x^'^'^'ii - . (24) 

We fix the s quark distribution to be 

/«^^(x)=0.3/«^°(x). (25) 
The gluons in the photon arc assumed to have hadron-like behaviour 

fa{x)^Al^^x^o^''{l-xfo^\ (26) 

As there are no data at a; ~ 1 we fix C^^^ — 1 and Cq^° = 3 as suggested 
by counting rules [28, 29] based on dimensional arguments. Thus we are left 
with six free parameters. 

Other, more flexible, forms of the x dependences at the starting scale 
were also investigated. No substantial improvement in the description of the 
data was observed, however the errors on the fitted parameters as well as 
the correlations were increased. Therefore we chose to present here results 
obtained with a minimal number of free parameters, as described above. 



4 Data sets 

For fitting the parameters we used 164 points of F2 measurements coming 
from e~^e~ reactions, 122 proton structure function data points from ep in- 
teractions and 24 points of dijet photoproduction. 
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4.1 data 

We have used all published data on the photon structure function F^, coming 
from LEP, PETRA, and TRISTAN^. 

The following LEP data have been included: 

• OPAL measurements from LEPl [32] and from LEP2 [33]. The LEPl 
data are in the range 1.86 < < 135 GeV^ (32 points), while the 
LEP2 data are in the range 9 < < 739 GeV^ (31 points); 

• L3 measurements [34] in the kinematic range of 1.9 < < 120 GeV^ 
(28 points); 

• DELPHI data [35] for = 12 GeV^ (4 points); 

• ALEPH measurements from LEPl [36] for = 9 9^ 20.7 and 284 
GeV^ (11 points), and from LEP2 [37] for = 17.3 and 67.2 GeV^ 
(16 points). 

The PETRA data are from PLUTO [38] at = 2.4, 4.3, 9.2 and 45.0 
GeV2, from JADE [39] at = 24 and 100 GeV^, and from TASSO [40] at 
Q2 = 23 GeV2. 

Finally AMY data [41] at = 5.8, 73.0 and 390.0 GeV^ and TOPAZ 
data [42] at = 5.1^ ig.Q and 80.0 GeV^ were used from TRISTAN. 

4.2 data 

As stated in section 12. 4^ we have used the Gribov factorization relation in 
order to produce indirect 'data' at low x from the proton structure func- 
tion data F2 measured by ZEUS [9]. The data have been moved to the 
values of the appropriate F^ data by using the ALLM97 parametriza- 
tion [43]. Only data with x < 0.01 and < 100 GeV^ were used. An 
additional systematic error of 3% was added to the measurement errors, to 
account for the uncertainty in the numerical coefficient used in the derivation 
of the indirect F2 'data'. The statistical and systematic errors were added 
in quadrature. 

4.3 Dijet photoproduction data 

The dijet photoproduction measurements were taken from the ZEUS exper- 
iment [11]. The cross section da/dx"''^ is measured in four bins of transverse 

■^The TPC27 data from PEP [30] have been excluded, as these data are considered to 
be inconsistent with other measurements [31]. 
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energy Et- 14-17, 17-25, 25-35, and 35-90 GeV. The jets are identified in the 
data by the clustering algorithm and is the transverse energy of the 
highest Et jet. 

5 Results 

The results of fits are presented for the following configurations: 

• F2 data, including the 'indirect' data from measurements, 

• F2 and dijet cross sections with the CTEQ5M parametrization of the 
proton PDFs [44], 

• F2 and dijet cross sections with the ZEUS-TR parametrization of the 
proton PDFs [45]. 

The masses of the heavy quarks used in the fit are rUc = 1.5 GeV,mfe = 
4.5 GeV and rrit = 174 GeV. 

The parameters are fitted to the data using MINUIT (release 96.03) [46] 
with MIGRAD and HESSE algorithms for error calculations. 

The results are summarized in Tabled Also added in the table are the 
values of the P^r degree of freedom and of the integral 1)211) representing 
the FG sum rule. 

Table 1: Parameters of the initial distributions at = 2 GeV^, per degree 
of freedom of the respective fits and the FG sum rule value at = 4 GeV^. 





no jets 


CTEQ5M 


ZEUS-TR 




4.04 ±0.26 


3.89 ±0.24 


4.45 ±0.29 


^PL 


1.22 ±0.20 


1.11 ±0.16 


1.90 ±0.23 


^HAD 


0.0645 ±0.0029 


0.0656 ±0.0029 


0.0647 ±0.0032 


j^HAD 


-1.17 ±0.006 


-1.16 ±0.006 


-1.16 ±0.006 


^HAD 


0.0173 ±0.0056 


0.0159 ±0.0049 


0.0271 ±0.0072 


rHAD 


-1.64 ±0.05 


-1.65 ±0.05 


-1.57 ±0.05 


XVNDF 


1.06 


1.53 


1.63 


FG SR 


1.13 


1.14 


1.04 



A general observation is that the addition of the dijet data has a minor 
effect on the overall fit results and their errors. There is though a noticeable 
deterioration in the value. Also the differences due to the usage of two 
different proton PDFs are relatively minor. Our preferred parametrization 
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is the one presented in the third column of Table ^ as this is the one which 
fulfils best the FG sum rule. We call it the SAL parametrization and it is 
presented in the following in more details. 



5.1 Comparison with F2 data 



a 



0.5 



: Q^=1.88GeV^ 
; q',p e [1 .86,1.9] GeV' 


= 2.55 GeV^ 
Q^pG [2.4,2.7] GeV ' 


= 3.63 GeV- 
^Q^ps [3.5,3.76] GeV ' 


1 Q^ = 4.95GeV^ 

[4.3,5.7] GeV ' 


■ = 7.43 GeV^ 
; q'^p g [6.5,8.5] GeV ' 


• = 10.3 GeV- 
Qexp^ [8-9.12] GeV ^ 


= 16.2 GeV- 
k Q^^pG [14.5,18] GeV ' 


' = 22.3 GeV- 
\ q'^p g [20.7,24] GeV ' 











10"' 10* 10"' 10"' 10"' 1 10' 10"' 10"' 10"' 1 10'" 10'" 10"' 10"' 1 10' 10"' 10"' 10"' 1 

X X X X 



; q' = 30.7 GeV' 
: qLp e [27,35] GeV ' 


- = 42.3 GeV^ 
_ q',p g [39.7,45] GeV ' 


- Q^ = 67.1GeV- 
!^Q^^^G^^76^ 


- Q^ = 89.4GeV- 
AQexp <^ [80,100] GeV ' 










: Q^=127GeV^ 
■ q',p g [120,135] GeV ' 


Q- = 284 GeV^ 


Q' = 390 GeV' 


\ Q- = 780GeV^ 
-\ Qf,p = 780GeV^ 











b1.5 
1 

0.5 



_ 2 
1 




10 10 10 10 1 10 10 10 1 10 10 10 1 10 10 10 1 

X X X X 

Figure 2: The SAL expectation for Q^) as a function of x at selected 

Q"^ values as denoted in the figure. The plotted data (dots for F2 measured 
directly and triangles for F2 deduced from Ff) are from the range Ql^^ 
presented in the figure. Note that the x range in the upper and lower plots 
is different. 



In Figure |21 we compare the F2 obtained with the SAL parametrization 
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together with the data, as a function of x in bins of Q^. The real data 
and the ones deduced from F2 are shown with different symbols. To limit 
the number of plots without loss of information, the data are shown within a 
range of Q^, while the corresponding curve is calculated for the average of 
that range. The shaded error band is calculated according to the final error 
matrix of the fitted parameters as returned by MINUIT. The uncertainty 
becomes smaller with increasing Q^, due to the expected loss of sensitivity 
to the initial conditions of the evolution. 

The inclusion of the F2 data in the fit strongly constrains the uncertainty 
on the gluon distribution. As an exercise, we increased the systematic un- 
certainty on the F2 'data' generated from from 3% to 10% and repeated 
the full fit. Small changes in the hadronic component of the photon were 
observed, leading to barely noticeable change in the PDFs in the range of x 
and investigated here. The most significant change was in the uncertainty 
on the gluon content at low and lowest x, were the error band increased 
by factor two. 




1 10 10' ,0 10 10% , 

Q [GeV ] Q [GeV ] 



Figure 3: dependence of for the x ranges shown in the plots. Black 
dots are F^ data and the line is F2 averaged over the depicted x range. 

The dependence of F2 is shown in Fig. 121 for different x ranges. This 
time only the F^ data are shown. Again each plot contains data from a 
range of x values, while the curve corresponding to the SAL parametrization 
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is averaged over the x range in the following way: 



X2 



X2 — Xi 



dxF^{x,Q^) . 



(27) 



A good agreement between data and the results of the fit is observed. 
5.2 Comparison with dijet data 

The dijet photoproduction cross sections are compared to calculations ob- 
tained using the Frixione-Ridolfi code with jets identified by the fcx-clustering 
algorithm, the same as used by the ZEUS collaboration to analyse the data. 
The predictions for the SAL photon parametrization and the ZEUS-TR pro- 
ton PDFs are presented in Fig.|31 Also shown are the individual contributions 
of partonic cross sections induced by gluons in the photon and gluons in 
the proton Gp. 




Figure 4: Dijet cross-section for ZEUS-TR proton x SAL photon — full line. 
Other contributions: Gp x (7 + q^) — shaded; (g^ + Gp) x G^ — dotted line. 
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With the exception of the lowest bin, the calculations tend to under- 
estimate the measured cross sections, especially at x < 0.4. It should be 
noted, however, that within the uncertainty^ of the NLO calculation there is 
no clear evidence for disagreement. Nevertheless, in this study we aimed at 
improving the agreement by adjusting the gluon content of the photon. How- 
ever, only in the range 0.1 < x < 0.2 the gluons from the photon contribute 
roughly half the cross section. In all other bins the contribution of the gluons 
in the proton dominates. The dijet data are therefore more sensitive to the 
gluons in the proton than to the gluons in the photon. It would therefore be 
tempting to use these data in the global fits of parton densities of the proton. 
However the present uncertainty on the NLO calculations of the dijet cross 
section is too large to reliably extract the gluon density in the proton. 

5.3 Parton distribution functions 

The SAL parton distribution functions in the photon are shown in Fig. 
The features to be noted are the behaviour of quarks at large x, typical of 
the point-like contribution of the photon, and the dominance of the gluon 
distribution at low x. 

A comparison between the SAL PDFs and the other available NLO DIS7 
photon parametrizations, GRV [3], GRS^ [6], and CJK [7], is shown in Fig. IHl 
for = 2.5 GeV^ and in Fig. [7| for = 100 GeVl 

At low there are big differences between the various PDFs^. They 
are especially pronounced for x < 10~^, where no F2 data are available 
and the result is subject to additional theoretical assumptions. The SAL 
parametrization has the lowest gluon distribution down to x ~ 10~^, below 
which value we observe a steep rise, steeper than in the other PDFs. At high 
Q^, where the sensitivity to initial conditions is diminished, there are still 
noticeable differences. 

To further compare the various parametrizations, the momentum fraction 
carried by the gluons and the light quarks are summarized in Table El They 
can be used to assess the fulfilment of the FG sum rule. According to the 
FG sum rule, the value of the sum G+S' should be ^ 0.9 at Q'^ = 2.5 GeV^ 
and ^ 1 at = 4 GeV^. The SAL parametrization is found to be closest 
to these values. 

^Changing the factorization/renormalization scale in the range between 0.5E^ and 2E^ 
may change the cross sections by about 20% [11]. 

^This parametrization uses Fixed Flavour Number Scheme (FFNS), where only u,d 
and s PDFs exist. 

non-vanishing 6-quark density at = 2.5 GeV^ is a feature of the CJK 
parametrization. 
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Table 2: Parton momentum fractions at = 2.5 GeV^ and 4 GeV^. E' 
= 2{p^+Pd + Ps) for the NLO PDFs. 





g2 = 2.5 GeV^ 




Q' = 


4GeV^ 






GRV 


GRS 


CJK 


SAL 


GRV 


GRS 


CJK 


SAL 


G 


0.47 


0.29 


0.87 


0.06 


0.50 


0.32 


0.90 


0.12 


d 


0.13 


0.07 


0.14 


0.10 


0.13 


0.08 


0.15 


0.10 


u 


0.25 


0.24 


0.21 


0.26 


0.28 


0.27 


0.24 


0.28 


s 


0.06 


0.04 


0.04 


0.07 


0.06 


0.05 


0.05 


0.07 


S' 


0.88 


0.70 


0.77 


0.85 


0.96 


0.78 


0.86 


0.92 


G+S' 


1.34 


0.99 


1.64 


0.91 


1.46 


1.10 


1.76 


1.04 
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6 Conclusions 

The F2 measurements supplemented by the low x expectations for F2 based 
on the measurements of F2 and Gribov factorization, and by dijet photopro- 
duction data have been used to extract a new NLO parametrization of the 
parton distributions in the photon. 

A good description of the data is obtained with the new parametrization, 
with the exception of the dijet data which turned out to be mostly sensitive 
to the gluon distribution in the proton. The obtained parton distributions 
fulfil the Frankfurt-Gurvich sum rule which was not imposed in the fit. 
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Figure 7: Comparison of SAL to other NLO parametrizations at 
100 GeVl 



Acknowledgements 

This work was partially supported by the Israel Science Foundation (ISF). 
One of us ("W.S.) greatly acknowledges discussions with Michal Praszalowicz 
and Krzysztof Golec-Biernat. 



References 

[1] See e.g. R. Devenish and A. Cooper-Sarkar, Deep Inelastic Scattering, 
Oxford University Press (2004). 



18 



[2] See e.g. M. Krawczyk, A. Zembrzuski and M. Staszel, Phys. Rep. 345 
(2001) 265; A. De Roeck, Eur. Phys. J. C33 (2004) S394. 

[3] M. Gluck, E. Reya, A. Vogt, Phys. Rev. D45 (1992) 3986, Phys. Rev. 
D46 (1992) 1973. 

[4] R Aurenche, J.-Ph. Guillet, M. Fontannaz, Zeit. Phys. C64 (1994) 621. 

[5] L.E. Gordon and J.K. Storrow, Nucl. Phys. B 489 (1997) 405. 

[6] M. Gluck, E. Reya and I. Schienbein, Phys. Rev. D 60 (1999) 054019; 
Erratum-ibid. D 62 (2000) 019902. 

[7] F. Cornet, P. Jankowski and M. Krawczyk, Phys. Rev. D 70 (2004) 
093004. 

[8] A. Levy, Phys. Lett. B404 (1997) 369. 

[9] ZEUS Collaboration, J. Breitweg et al, Eur. Phys. J. C7 (1999) 609; 
Eur. Phys. J. C21 (2001) 443. 

[10] V.N. Gribov, L.Ya. Pomeranchuk, Phys. Rev. Lett. 8 (1962) 343. 

[11] ZEUS Collaboration, S. Chekanov et al, Eur. Phys. J. C23 (2002) 615. 

[12] A. Donnachie and P. Landshoff, Phys. Lett. B296 (1992) 227. 

[13] W. Furmanski, R. Petronzio, Zeit. Phys. Cll (1982) 293. 

[14] S. Frixione, Z. Kunszt and A. Signer, Nucl. Phys. B467 (1996) 399; 
S. Frixione, Nucl. Phys. B507 (1997) 295. 

[15] E. Witten, Nucl. Phys. B120 (1977) 189. 

[16] W.A. Bardccn and A.J. Buras, Phys. Rev. D20 (1979) 166[Erratum- 
Phys. Rev. D21 (1980) 2041]. 

[17] R.J. DeWitt, L. M. Jones, J.D. Sullivan, D.E. Willen and H.W. Wyld, 
Phys. Rev. D19 (1979) 2046[Erratum- P/iys. Rev. D20 (1979) 1751]. 

[18] G. Curci, W. Furmanski, R. Petronzio, Nucl. Phys. B175 (1980) 27; 
W. Furmanski, R. Petronzio, Phys. Lett. 97B (1980) 437. 

[19] M.A.G. Aivazis, J.C. Collins, F.l. Olness, and W.-K. Tung, Phys. Rev. 
D50 (1994) 3102. 

[20] J.C. Collins, Phys. Rev. D58 (1998) 094002. 



19 



[21] R.S. Thorne and R.G. Roberts, Phys. Rev. D57 (1998) 6871. 

[22] W.-K. Tung, S. Kretzer and C. Schmidt, J. Phys. G. 28 (2002) 983. 

[23] V.M. Budnev, I.F. Ginzburg, G.V. Meledin and V.G. Serbo, Phys. Rep. 
15 (1974) 181. 

[24] For a recent discussion see R. Thome's talk at the DIS05 workshop, 
[hep-ph/05062^ 

[25] G.A. Schuler, T. Sjstrand, Zezt. Phys. C68 (1995) 607. 

[26] L.L. Frankfurt, E.G. Gurvich, "hep-ph/9 505"406l Phys. Lett. B386 

(1996) 379;J. Phys. G. 22 (1996) 903. 

[27] W. Slominski, to be pubhshed. 

[28] R. Blankenbecler and S.J. Brodsky, Phys. Rev. DIO (1974) 2973. 

[29] G.R. Farrar and D.R. Jackson, Phys. Rev. Lett. 35 (1975) 1416. 

[30] H. Aihara et ah, Phys. Rev. Lett. 58 (1987) 97; H. Aihara et al., Zeit. 
Phys. C34 (1987) 1. 

[31] S. Albino, M. Klasen and S. Sldner-Rembold, Phys. Rev. Lett. 89 (2002) 
122004. 

[32] OPAL Collaboration, Zeit. Phys. C61 (1994) 199; Phys. Lett. B412 

(1997) 412; Zeit. Phys. C74 (1997) 33; Eur. Phys. J. C18 (2000) 15. 

[33] OPAL Collaboration, P/i2/s. Lett. B411 (1997) 387; Eur. Phys. J. C18 
(2000) 15; Phys. Lett. B533 (2002) 207. 

[34] L3 Collaboration, P/iys. Lett. B436 (1998) 404; Phys. Lett. B447 (1999) 
147; Phys. Lett. B483 (2000) 373. 

[35] DELPHI Collaboration, Ze2t. Phys. C69 (1996) 223. 

[36] ALEPH Collaboration, Phys. Lett. B458 (1999) 152. 

[37] ALEPH Collaboration, Eur. Phys. J. C30 (2003) 145. 

[38] PLUTO Collaboration, Phys. Lett. B142 (1984) 111; Nucl. Phys. B281 
(1987) 365. 

[39] JADE Collaboration, Phys. Lett. B121 (1983) 203; Zeit. Phys. C24 
(1984) 231. 

20 



[40] TASSO Collaboration, Ze^^. Phys. C31 (1986) 527. 

[41] AMY Collaboration, Phys. Lett. B346 (1995) 208; Phys. Lett. B400 
(1997) 395. 

[42] TOPAZ Collaboration, F/iys. Lett. B322 (1994) 447. 

[43] H. Abramowicz and A.Levy, The ALLM parameterization of atotil * p) ■ 
An update., Preprint DESY-97-251 |hep-ph/9712415|, DESY, 1997. 

[44] CTEQ Collaboration, H.L. Lai et al, Eur. Phys. J. C 12 (2000) 375. 

[45] ZEUS Collaboration, S. Chekanov et al, Phys. Rev. D 67 (2003) 012007. 

[46] For the fits the MINUIT package is used: MINUIT (release 96.03): 
F. James and M. Koos, Computer Phys. Comm. 10 (1975) 343. 



21 



